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Abstract. The heating of the solar corona is one of the big questions in astrophysics. 
Rapid pulses called nanofiares are among the best candidate mechanisms. The analysis 
of the time variability of coronal X-ray emission is potentially a very useful tool to 
detect impulsive events. We analyze the small-scale variability of a solar active region 
in a high cadence HinodefKRT observation. The dataset allows us to detect very small 
deviations of emission fluctuations from the distribution expected for a constant rate. 
We discuss the deviations in the light of the pulsed-heating scenario. 



1. Introduction 

The investigation of the heating mechanisms of the confined coronal plasma is still 
under intense debate. It is widely believed that the energy source for coronal heating 
is the magnetic energy stored in the solar corona. An unsolved problem is how this 
magne tic energy is converted into thermal energy of the confined coronal plasma. As 
Parker proposed in 1988, rapid pulses called nanofiares are among the best candidate 
mechanisms of magnetic energy release. Nowadays a challenging problem is to obtain 
evidence that such nanofiares are really at work. If small energy discharges (nanofiares) 
contribute in some way to coronal heating, they could be too small and frequent to be 
resolved as independent events. In this case, we would need to search for indirect 
evidence. 

The idea of this work is that, if the solar corona emission is sustained by repeated 
nanofiares, locally the X-ray emission may not be entirely constant but may show 
variations around the mean i ntensity. So nanofiares m a y leave their signature on the 
light curves. Many authors dShimizu & Tsunetalll997l: Vekstein & Katsukawal 120001 : 



Katsuka wa & Tsu neta 2001 ; Katsu kawall2003l : ISakamoto et al . 2008) pointed out that a 
detailed analysis of intensity fluctuations of the coronal X-ray emission could give us 
information on these smallest flares. Following this hint we use this approach for the 
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first time on Hinode data, searching, with statistical analysis, for small but systematic 
variability in noisy background light curves and their link to coronal heating models. 



2. Instruments and Data 



The Japanese Hinode satellite was launched in 2006 (Ko sugi et alJl2007T) . One of the 
instrumen ts it carries is the X-Ray Telescope (X RT), a high resolution grazing incidence 
telescope (iGolub et alOool : iKano et a l.l l2008f) that images the Sun through nine X-ray 
filters, using two filter wheels. 

In our work we analyze the active region (AR) 10923, observed when it was lo- 
cated close to the disk center on November 14, 2006. The data set consists of 303 
images of 256 x 256 pixels, taken in the ALpoly filter with an average cadence of about 
6 s, for a total coverage of the observation of about 26 minutes. 



3. Data Analysis 

In this work we analyze the intensity fluctuations of a high cadence XRT observation, 
searching for signal on top of statistical fluctuations, and its link to coronal heating 
model. Figure |2]is a flow chart which describes the concept of our analysis. Let us con- 
sider the light curve of a pixel in the AR and let us build the histogram of the amplitude 
of the fluctuations with respect to the average emission trend. If the count rate is strictly 
constant, the histogram of fluctuations is trivial. However, fluctuations are intrinsic to 
the finite photon counting. If we add the expected photon noise and the count rate is 
high, we expect an almost Gaussian distribution of the fluctuations, symmetric around 
zero, i.e., we expect an equal amount of positive and negative fluctuations. We expect 
the same distribution if the count rate is linearly and moderately changing instead of 
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Figure 2. Schematic concept of this work. Consider the light curve of a valid 
pixel. If the count rate is constant, the histogram of fluctuations is trivial. If we add 
the expected photon noise and the count rate is high we expect a broad Gaussian 
distribution of the fluctuations. We expect the same distribution if the count rate is 
linearly changing. 



being constant. Even for a large curve with a constant rate we expect a broad distri- 
bution of fluctuations due to photon noise. If the count rate is high the distribution is 
approximately Gaussian. We expect exactly the same distributions also if the count rate 
is linearly changing instead of being constant, as sketched in Figure |2] 

Along this line, for each pixel we create the distribution of intensity fluctuations 
of the light curve. We normalize the amplitude of the fluctuations to the photon noise 
cr p , expected from determining the original photons impacting on the detector from the 
measured readout. The conversion from DN to photon counts is through the conversion 

(2) 

factor K for the z-th filter-band: 

cr v = Jk?\T)I , (1) 

where T is the electron temperature and Io is the measured DN. 
The conversion factor is obtained as: 

(2) J [(hc/A)/(51 x 3.65eV)] P(A, T) m{A) dA 

K: = r , (2) 

fP(A,T) m (A)dA 



where P(A, T) is the emissivity as a function of wavelength A and electron temperature 
T, and rg(A) is the telescope effective area. For the mome nt, we assume a temperature 
log T = 6.3 dKano & Tsunetall 19951 : IRarukage et al.ll201lh . 
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Figure 3. Top: Distribution of the fluctuations of a pixel light curve, with re- 
spect to the linear fit to the light curve (solid line) in comparison with a symmetric 
Gaussian-like distribution (dashed line). Bottom: Residuals of the fluctuations with 
respect to the Gaussian curve. 



Since we are interested in low amplitude systematic variations, we first remove 
pixels with low signals and those showing all kinds of high or slow amplitude system- 
atic variations, such as: 

• Pixels with spikes due to cosmic rays 

• Pixels that show micro-flares or other transient brightenings 

• Pixels that show slow variations due, for instance, to local loop drift or loop 
motion. 

At the end of this data screening we are left with about 56% of total pixels. The 
light curves of all these pixels can be fit with linear regression. An an example of our 
analysis, in Fig.[3]we show the result for one pixel (1,1). 

3.1. Results 

For each pixel we build the histogram of the X-ray intensity fluctuations with respect 
to the linear fit line. Figure [3] shows a sample distribution for one pixel compared 
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Modeling: small flare-like perturbations 



i 



I 




t 



t 





Median < 



Median < 



Figure 4. Schematic concept of expectations from light curves perturbed by trains 
of small flare-like events. Although the mean emission would not change, the distri- 
bution of the fluctuations would become asymmetric, with a negative median even 
in the presence of significant photon noise. 

to a symmetric Gaussian distribution. By inspecting the figure, we realize that there 
is a slight excess of the number of negative fluctuations over the number of positive 
ones. This excess makes the histogram slightly asymmetric toward the left side of the 
distribution. The excess is present but not very significant on most of the single pixel 
signal over the noise (as we can see in Figure [3]). 

From preliminary analysis, the excess becomes more significant when we apply 
the analysis to selected subregions localized in different part of the active region and 
even more when we sum over the whole region. 



4. Conclusions 

In summary, histograms of light curve fluctuations measured at high cadence show 
asymmetric distributions with median lower than zero. Asymmetries seem to appear at 
all scales from single pixels, sub-regions, to the whole active region at high confidence 
level. We propose that the asymmetry of the distributions of the fluctuations could be 
compatible with random se quences of s mall fast rise and slow decay pulses. This is typ- 
ical of flare-like evolution (Realel l2010h . Figure [4] sketches the concept of this scenario: 
the presence of trains of small flare-like events would make the distributions of fluctu- 
ations asymmetric in the observed direction. The median would be negative because 
the light curve will spend more time below than above the average trend, i.e., the mean 
for a constant light curve. Preliminary modeling using Monte Carlo simulations ap- 
pears to support this scenario. In our opinion, if confirmed, these results would support 
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a widespread storm-of-nanoflares scenario. We emphasize that this analysis requires 
high-cadence XRT observations. Improved diagnostics would be obtained if at least a 
few images were available in different filter bands. 
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